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Abstract 

The spectra of charged particles in the energy variable !^ = log Ij x ("the hump 
backed plateau") are analyzed in terms of moments in order to investigate the 
sensitivity to the running of in particular towards small scales, and to further 
test thereby the local parton hadron duality picture. The modified leading log 
approximation (MLLA) of QCD provides a very satisfying description of the 
e"'"e~ data in the full range of cm. energies from about 1 up to 91 GeV. 
The same calculation but with fixed coupling yields moments which are 
inconsistent with the data already close to threshold. In particular the soft 
particles {E <1 GeV) refiect directly the running of a^. 
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1. Introduction 

The predictions from perturbative QCD on the intrinsic structure of jets based on the idea 
of a soft confinement mechanism [1, 2] have been remarkably successfuL In particular, the 
prediction of the "hump backed plateau", i.e., an approximately Gaussian distribution 
in log(Pjet/_Phadron) with a supprcssiou of low momentum particles from the soft gluon 
interference [3, 4, 5, 6] has been confirmed experimentally. This demonstrated the close 
similarity of parton and hadron spectra, provided the parton cascade is evolved down 
to a cutoff scale of the order of the hadronic masses Qq ~ nih- In this way, an almost 
quantitative description of the momentum spectra in the PETRA/PEP energy range (see, 
e.g., [7]) and at LEP[8, 9, 10, 11] has been obtained ("Local Parton Hadron Duality" [12]). 
The first data from IIERA[13, 14] confirm this picture as well. 

In this paper we investigate the sensitivity of the energy spectrum of final state particles 
to the running of the coupling constant a^. The effect of the running is strongest for 
soft particles and at low energies and its verification can therefore provide an additional 
hint on the applicability of perturbative QCD towards small scales. To this end, we use 
a consistent kinematical scheme to relate partons and hadrons. We analyze the spectra 
in terms of moments for which explicit analytical expressions can be derived for fixed and 
running coupling and we extend the analysis down to low energy of order 1 GeV. 

An appropriate theoretical framework is the evolution equation for parton densities in 
the modified leading log approximation (MLLA) (for a recent review see [2]). It determines 
the parton densities and moments in absolute magnitude by the boundary condition at 
threshold in terms of two parameters, the QCD scale A and the pt cutoff Qq. This is 
different from deep inelastic scattering, where an unknown parton distribution has to be 
provided at a particular scale. The asymptotic predictions at high energies agree with 
those from the double log approximation (DLA) which neglects the energy-momentum 
constraints[3, 4]. The Next-to-Leading corrections of relative order y^aj for moments are 
found large[5]. The MLLA results[6] sum all orders of y^aj which allows the extrapolation 
to threshold but not all contributions beyond Next-to-Leading order are taken into account. 
More details of this study will be published elsewhere, see also [15]. 



2. The theoretical framework 



MLLA Evolution Equations. The evolution equation of single parton inclusive energy 
distribution in the framework of MLLA is given by: 



d 



-xL'|(x,logP0) 



E 



Ztt 



d log r< - 

o C=q,q,g 

with the boundary condition at threshold P0 = Qq: 

xD^{x^ log Qo) = 5(1 — x)5 
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Here P and denote the primary parton energy and jet opening angle respectively, x is 
the energy fraction carried by the produced parton, $j (z) are the parton splitting kernels 
and A, P, C label quarks, antiquarks and gluons. The QCD running coupling is 
given by its one-loop expression as{kt) = 2tt /blog[kt/ A) with b = (llA^c — 2n/)/3, A the 
QCD-scale and Nc and nj the number of colors and of flavors respectively. The scale of 



2 



the coupling is given by the transverse momentum kf ~ z[l — z)PQ. The shower evolution 
is cut off by Qo, such that kf > Qq. 

The integral equation can be solved by Mellin transform: 

D^Y, A) = /' —x'^[xD{x, ¥)] = r d^e-^'^D{^, Y, A) (3) 
Jo X Jo 

with Y = log ~ log A = log = log and parton energy k. 

In flavor space the valence quark and (it) mixtures of sea quarks and gluons evolve 
independently with different "eigenfrequencies" . At high energies, the dominant con- 
tribution to the inclusive spectrum comes from the "plus"-term, which we denote by 
-Dc^(y, A) = _D^(y, A). In an approximation where only the leading singularity in cu-space 
plus a constant term is kept, it satisfies the following evolution equation: 

where a = llA^c/3 + 2nf /3N^. Taking a = and dropping the r.h.s would yield the DLA 
evolution equation. The analysis of this paper is based on Eq. (4). 
By defining now the anomalous dimension according to: 

D^Y, A) = D^{Yo, A) exp (^J^ dy-f^[as{y + A)]^ (5) 

the evolution equation for the inclusive spectrum can be written as a differential equation 
for as follows: 

+ 7c^)7c^ - ^^^^ = -f3{a,)-^-f^ - 2r]{(jj + -f^) (6) 
ZTT dag 

where f3{as) = -^as{Y) ~ — and r] = a^^/SNc = aas/Air. In general this nonlinear 
equation has two roots and the solution consists of a superposition of two terms of the 
type (5). 

Moment analysis for running as- Explicit expressions with parameters Qq and A can be 
derived for the multiplicity A/" and the normalized moments < (^'(K, A) >= / d^^'^ D[^^Y^ X) / J\f 
or the cumulant moments Kq[Y^ 6]; they are related by Hi =< ^ >= K2 = cr^ =< 
— 1^)^ >, K3 =< ((^ — >, K4 =< ((^ — > — 3(t'*, . . . ; also one introduces the reduced 
cumulants kg = Hqja^, in particular the skewness s = k^ and the kurtosis k = k4. The 
cumulant s be obtained from: 



^^,{Y,X) =(-^y log DUY,X) , (7) 
also Me = -De^l^-o- At high energies, one term of the form (5) dominates and one obtains 

(8) 



Kg{Y) = Kg{Yo) + 7c^[«s(y)] 



uj=0 
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This equation shows the direct dependence of the cumulants on as{Y). For fixed a^, 
for example, one obtains directly Kq[Y) oc Y for high energies. In our applications, we 
take into account the boundary conditions at threshold 1^ = 0, eq. (2), which yields 
K,{Yo = 0) = 0. 

The eqs. (5) and (8) are well suited to discuss the dependence on the number of 
fiavors nj. The simplest procedure would consist in taking 7t^(y) under the integral as 
step function with nj as the number of open fiavors at energy y. In leading order of 
\/Y (, i.e., in the DLA), nj enters only through the running coupling as{y^nj)^ at non 
leading level also through the parameter a (from the splitting function $5). Restricting 
our discussion to the leading term we note that in the original eq. (1) the scale of is 
kf ~ zil — z)P < P/A and the maximal scale P appearing in eqs. (5) and (8) originates 
from a kinematical approximation. In our comparison with data we choose either nj = 3 
throughout as in previous analyses or we include the effect of heavy quarks in eq. (5) and 
(8) with nj changing for momenta Pq = AniQ. A more accurate treatment of scales would 
have to take into account two-loop contributions besides the above kinematical effects and 
possibly explicit mass dependence of the coupling[16]. 

The evolution equation for D^^[Y^X) can be analytically solved and one obtains[2] a 
linear superposition of two hypergeometric distributions. Moments of order q can then be 
obtained through the formula < > = [( — -^^^(^5 -^)] L=o/'^ written in 

general as: 

<e>=;i^E (^)(^i4'^ + ^2i??^) (9) 

where iVi, iVa, L^k^ and i^^"^ are known functions of a, 6, K + A and A. So they depend on 
the two parameters Qq and A[6]. We found, as discussed below, the best agreement for 
the "limiting spectrum", where the two parameters coincide, Qq = A, i.e., A = 0. In this 
case, the formulae simplify and the moments can be expressed[6] in terms of the parameter 
B = a/b and the variable z = ^JlQNcY/b as: 

= pI,'\b + 1,5 + 2,.)+ 2 ^£±iMp(^)(i^ + 1,5 + 2,.) (10) 
where Pq'^ and P^''^ are polynomials of order 2{q — 1) in z. These moments have an 



expansion < > /K* = Cq'^ + c^'^/V^ + • • • ^- The higher orders of this series are 
still numerically sizeable at LEP energies (10% contribution from NNLO to ^ and cr^) 
and increase towards lower energies. We therefore included the full expression, also for 
P/'^; explicit results for the moments for g < 3 can be found in [6] and for g = 4 in [15]. 
The average multiplicity of partons (zeroth-order moment of the inclusive momentum 
spectrum) for the Limiting Spectrum is given by: 



MlS = r(P) y iB + ^iz) . (11) 

This expression does not fulfill the boundary condition M{0) = 1. Therefore in our fits we 
added a constant term"*. 



^We checked that ^ and agree in NLO and s in LO with the result obtained by Fong and Webber[5], 
i.e., agreement of < > in NLO. 

■*The general solution considered in [6] yields dD^ /dY\Y-Q= ~B/X; therefore the multiplicity for 
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Moment analysis for fixed as- In order to study the sensitivity of the inclusive energy 
distribution to the running of the couphng, let us build a new model based again upon the 
MLLA but with the coupling or, equivalently, the anomalous dimension 70 = ^JGo^Jtt, 
kept fixed as a free parameter replacing A. In this case, the evolution equation for the 
Mellin transform, D^^[Y^X)^ eq. (4), becomes: 



12) 



One can then perform explicitly the inverse Mellin transform 



d^dY 



dV 



Through the simple transformation: D[^^ A) = -D((^, A) exp[— 2?7(y — (^)], one gets rid 
of the linear terms and -D((^, A) obeys eq. (13) with = 0, i.e., the DLA equation. Using 
its known solution[2], we find for the inclusive energy spectrum in MLLA with fixed-a^ 



D{i, F, A) = 70^^/1 ( 270^/^(1^ - ) e-'^^"--^^ 



(14) 



The exact form of cumulants can also be obtained in a straightforward way from the 
evolution equation for the anomalous dimension 7^^ given in eq. (6). For fixed a^, the first 
term on the r.h.s. drops out; we obtain an algebraic equation with two solutions 7e^_|_ and 
we find: 



D^{Y,\) 



7c. ± 



7u 



7 



+ 



7 



7 



7a 



cu + 2r7)± V(^-2^) +470 



(15) 



(16) 



The results for the moments can be obtained easily from eq. (7); for example, we find for 
the average multiplicity Mjix'- 



AT, 



fix 



cosh (70^) + — sinh (70^) 
7o 



exp [—rjY) 



(17) 



with 7o = ^7q + 77^; it reduces for r] = to the well-known DLA result cosh(7oy"), and 
for the first moment 



fix 



[A cosh (70F) + B sinh (70F)] 



exp [—rjY) 



18) 



fix 



where A = Y{^q + ri^)/2^Q and B = {2^Q{riY — 1) — ri^)/2^Q. By choosing the + term 
in eq. (16), one selects the leading high energy contribution which also dominates the 
"limiting spectrum" for the running model. In this limit, according to eq. (8), one gets: 



AT^ 



(19) 



general A has an unphysical behavior near threshold (i7V(Y)/(iY|y_Q< with 7V(0) = 1. For A = 0, 
#(0) = and rfA'(Y)/rfY|y^o> 0. 
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3r] 1 , 3{4:r]^ - 7^ 

70 V^oY 

With ?7 = 0, we obtain back the DLA results 



fix — / _ -.^ 1 f^fix — 9_ -.^ (21j 

70 VToY 7o7o >^ 



3. Relating parton and hadron spectra 

In the theoretical calculations, we take the partons as massless, but introduce a pt cutoff 
Qo for regularization; on the other hand, the observable hadrons are massive. One can 
make the assumption that the cutoff Qq is related to the masses of hadrons. As a first 
stage in the discussion of charged particle spectra, one can relate Qq to an effective hadron 
mass[12], then for both partons and hadrons E > Qq. Furthermore we relate the spectra 

Eh — = Ep — ; Eh = Ep > Qo (22) 
dph dpp 



at the same energy E or = log P/E (not momentum!), where Eh = ^JpI ~'r Qo and 
Ep = Pp. With this choice both spectra vanish linearly at the same E oi ^e- To see this, 
we note that the hadronic data behave approximately like 

i?|^~e-^/^° or E^c:,8nQoiE-Qo)e-^^^° (23) 

for sufficiently small energies (see Figure 3a below). The same linear decrease with en- 
ergy towards the kinematic limit is also followed by the parton distribution, Edn/dp ~ 
logE/Qo ~ (£" — Qo)/Qo (eq. (14)) and corresponds to the limited phase space with 
Pp > Qo- 



4. Results 

Moments. The moments < ^'^ > are determined from the spectra Edn/dp vs. ^e after 
appropriate transformation of the measured Xp = 2p/yG spectra of charged particles using 

E = \/p^ -\- Qo ^iid therefore depend on Qo. For the unmeasured interval near ^e — Y 
(small momenta) a contribution was found by extrapolation. The error for the moments 
includes the errors on the central values of ^e in each bin, taken as half the bin-size. In 
the same way we also obtain the multiplicity JVe as integral over ^e of the full spectrum 
Edn/dp; its difference to the usual particle multiplicity h depends on the cm. energy, 
ranging from a decrease of 30% at = 3 GeV to a decrease of 10% at LFP energy^. First 
we perform a fit of the two parameters Qo and A by comparing the moments determined 
for a selected Qo with the theoretical predictions from eq. (9) of different A. The best 
agreement is obtained for 

^ A ~ 270 MeV (24) 

®the MARK I data point at ^/s = 4.03 GeV shows an anomalous decrease up to 50%, which may be 
partly related to charm thresholds effects. 
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and we estimate AQq ~ AA ~ 20 MeV. The cumulants up to order 4 together with 
the corresponding MLLA predictions are shown in Figure 1. We also analyzed in the 
same way results from AD0NE[17] at -^i = 1.6 GeV which however reports only charged 
pions and therefore cannot necessarily be considered on the same ground as the other 
data. Nevertheless it is interesting to observe that the results follow smoothly the trend 
of the higher energy data, in particular = 0.53 it 0.07 and cr^ = 0.07 it 0.02. Errors 
of the cumulants are taken as the sum of statistical and systematic errors. The latter 
particularly affects the overall normalization and then the average multiplicity; higher 
order cumulants are on the contrary independent of the overall normalization. The LEP 
average has been obtained by following the PDG procedure[18] for weighted averages. The 
MLLA predictions are calculated for the number of flavors nj = 3. Including the heavy 
flavors at scale P/4 as discussed above would lower the prediction of cr^ at LEP by two 
std. and by less than one std. at lower energies; likewise the effect on 5, and k is 
negligible. 

The result — A ("limiting spectrum") follows, as an increase of A = log((5o/A) 
would predict lower values for ^e but larger values for cr^. These changes cannot be cured 
by a lowering of Qo, which would shift both ^e and cr^ downwards. This can be seen from 
the first two terms in the expansion of Ce and (7^: ^ = Y/2 + ab{VYTX - VA)/24 and 

= ((y + A)3/2-A3/2)/66-6y/96. Eor the multiphcity we write A^i^ = cipMLS + C2 with 
arbitrary normalization parameters for JVls in eq. (11) (the factor 4/9 is for the quark jet 
and 2 is for the two hemispheres) and the leading particle contribution. Eitting the lowest 
and highest energy data points we find Ci = 1.43 and C2 = 1.02. 

We have not been able to obtain a satisfying description for the fixed-a^ case. We show 
in Eig. 1 one illustrative example for the same Qq = 270 MeV and 70 = 0.64 (with nj = 
3). With this choice for 70 one obtains a good fit for the multiplicity, where we proceed 
as above and we write JVe = CipAfjix + C2 using eq. (17). In this case Ci = 3.10 and 
C2 = -1.24. Also the asymptotic slope for ^e is well reproduced. An adjustement of the 
absolute normalization of the moments ^e at a particular energy Yq is possible if the PQCD 
evolution towards lower energies and the normalization at threshold are abandoned. 

The predictions from MLLA with running and normalization at threshold are re- 
markably successful considering the fact that there are only two parameters, actually 
coinciding, for the four moments. A deviation is visible for ^e at the lower energies; this 
may indicate the infiuence of the leading valence quark which is neglected by restricting 
to the D'^ spectrum. Otherwise the predictions are confirmed at an almost quantitative 
level down to the lowest energy where charged particle spectra are measured. Also the 
log JVe data show a clear downward bending predicted from the running if the lower 
energy data points are included. The distinction from the fixed case would however 
require still higher energies. Excluding the lowest two data points of we find an overall 
X^NDF ~ 1.8. 

On the other hand, the same model but with fixed coupling only fits the data on JVe 
and the slope of ^e but fails otherwise. A fixed regime can be excluded already close 
to threshold. Indeed, if we suppose that the coupling is fixed for a certain energy interval 
near threshold < Y < Yi and only runs for Y > Yi^ then ^ would be shifted towards 
smaller values because of the very different evolution near threshold; for example for Yi 
= 0.3 (-^i ~ 0.7 GeV) the shift would be ~ 0.25, considerably away from the data. 
Also, if we restrict the discussion to higher energies, one cannot fit the higher moments 
with the same parameter 70 used for JVe and ^e^ as can be better seen from the rescaled 
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cumulants. 

Rescaled cumulants. In addition to the standard moment analysis performed in the pre- 
vious subsection, we also consider the rescaled cumulants Hqj!^. These quantities become 
energy independent in case of fixed coupling (see eq. (8)) at high energies and therefore 
exhibit more directly the difference to the case of running coupling. In Figure 2 we show 
the experimental results on these ratios which have a considerable y-dependence and are 
well reproduced by MLLA with running in the available energy range. The Figure also 
shows the effect of including the heavy quarks at scale P/4 as discussed above. It would 
be interesting to continue this type of studies with data from jets of higher energies, where 
the MLLA predictions for these ratios continue to show a remarkable dependence on Y . 

Low-energy behavior: Boltzmann factor and running-as effect. Let us consider the 
Lorentz-invariant distribution Edn/d^p as a function of the particle energy E for small 
energies (up to 1 GeV). Fxperimental data from low energy experiments show an exponen- 
tial decrease as in eq. (23) with Eq of the order of 150 MeV[19, 17, 20]. This behavior is 
often related to the Boltzmann factor in the thermodynamical description of multiparticle 
production as, for example, in Hagedorn's model[21]. 

Figure 3a shows experimental data on the invariant distribution at three different cm. 
energies. All three distributions tend to a common behavior at very low energies but 
deviations from a simple exponential are visible at larger cm. energies. In Figure 3b, 
we also show the theoretical predictions for the MLLA with running and fixed at 
= 3 GeV. The running-ag prediction has been obtained by using the distorted Gaussian 
parametrization as in [5] with the first four moments computed for the "limiting spectrum" 
with the same value of Qq. The fixed-a^ curve is obtained from the exact solution for the 
energy spectrum given in eq. (14) again with 70 = 0.64. Both theoretical predictions are 
then normalized to the experimental average multiplicity. 

It is remarkable that the MLLA with running is in agreement with experimental 
data even at such a low energy; the exponential decrease is well reproduced by the running- 
as effect, which mimics a thermodynamical behavior. On the contrary, the fixed-a^ model 
shows a fiatter behavior, inconsistent with experimental behavior. As for decreasing par- 
ticle energy E also the typical particle px is necessarily decreasing, the coupling as{pT/^) 
is rising in the running case. This yields the steeper slope in comparison to the fixed 
as case. Therefore the running of the QCD coupling is visible even in the energy range of 
few hundreds MeV in the particle spectra. 

The behavior near the kinematical limit E ^ Qq can be easily understood analytically 
in the DLA. The first iteration of the evolution equation of the energy spectra[2] (the 
Born term) yields a finite contribution for E ^ Qq and does not depend on cm. energy 
whereas the second term provides the rise of the spectrum for large E with increasing 
for both fixed and running as[15] in qualitative agreement with Figure 3a. For the ratios 
of invariant cross sections, I{E) = Edn/d^p^ we find for small E: 

IrunjE) ^ J2_ / 1 log E/Qo 

IjUE) ^hf^ V 2 A 

There is a large relative enhancement of this ratio from the running of at small E > Qq 
followed by a steep decrease for typical parameters (I2/67Q ~ 4.3 for 70 = 0.64, A < 1), 
as the pt integral of the Born term is dominated by small values of pt < E and therefore 
large coupling as{pt)- 
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5. Conclusions and final remarks 

The dependence on cm. energy of the experimental first four order moments of the 
inclusive log 1 j x distribution is well reproduced by the MLLA of perturbative QCD applied 
to the full cm. energy region down to threshold. We found that the best agreement with 
data is for the limiting spectrum, where Qq = A; the best value of Qq has been estimated 
to be around 270 MeV and corresponds to an effective mass of a charged particle. The 
zero-th order moment, i.e., the average multiplicity, requires two more parameters which 
fix the overall normalization for particle production and the leading particle contribution. 
The effect of the running is most pronounced near threshold where the variation of the 
coupling is strongest; it is also directly seen in the steepness of the slope of the invariant 
energy spectrum at low particle energies E < 1 GeV. Predictions of the same model 
but with fixed coupling are inconsistent with the experimental behavior in the full energy 
range. We take the evidence for the running effects at low energies E and y/s as further 
support for the local parton hadron duality hypothesis. 

The distribution of particles at low energies ^/s < 3 GeV is not just according to phase 
space: the running enhances configurations with coUimated partons. At the hadronic 
level such configurations correspond to hadronic resonances which are known to dominate 
the final states already at these low energies. It is therefore suggested to consider a duality 
between partonic and hadronic distributions already at low energies. The running with 
strong coupling at small pT seems to be related to the production of resonances. 

A word of caution is appropriate nevertheless. The MLLA results used here are based 
on certain high energy approximations which might not be fully justified at low energies 
(restriction to D'^ contribution, neglect of Next-to-Next-to-Leading effects). A deviation 
in at 3 GeV could hint towards a need of further contributions. Also two loop effects, 
whose importance increases towards lower energies, have been neglected in the calculation. 
At the price of introducing new parameters one could of course restrict the comparison 
with data to an energy range V > 1^ > above threshold. At any rate, the high 
energy approximations work well for energies > 3 GeV, where data on charged particle 
distributions have been presented. 

The higher moments q > 2 show a remarkable energy dependence in the energy range 
above the LEP experiments. This suggests that the analysis of hadronic jets along these 
lines at the TEVATRON at ten times higher energies or at future colliders should be able 
to test the proper running of towards much higher energies. 
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Figure Captions 



Fig. 1: The average multiplicity Me and the first four cumulants of charged particles' 
energy spectra Edn/dp vs. i-e., the average value ^e^ the dispersion cr^, the skewness s 
and the kurtosis A;, are shown as a function of Y = log(Y^/2(5o) for Qo = 270 MeV. Data 
points from MARK I[20] at = 3, 4.03, 7.4 GeV, TASS0[7] at = 14, 22, 35, 44 GeV 
and LEP at = 91.2 GeV (weighted averages of ALEPH[10], DELPHI[11], L3[9] and 
0PAL[8]). Predictions of the "limiting spectrum" (i.e. Qq = A) of MLLA with running 
as and of MLLA with fixed are also shown (for nj = 3). Predictions of the average 
multiplicity refer to the two-parameter formula JVe = Ci^2M + C2. 

Fig. 2: Rescaled cumulants Hqj as a function o{Y = \og{^/s|2Qo) with the correspond- 
ing predictions of the "limiting spectrum" of MLLA with running either without or 
with heavy fiavors included. Data as in Eigure 1. Eor fixed these quantities approach 
constant values at high energies. 

Fig. 3: Invariant distribution Edn/d^p as a function of the particle energy E for Qq = 270 
MeV. a): experimental results from MARK I[20], TASS0[7] and 0PAL[8]; b): comparison 
of experimental data from MARK I[20] with predictions from MLLA with running and 
with fixed at the same energy. 
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